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ABSTRACT 



The theory and equations for forecasting the wind drift of sea ice, 
as previously developed by several researchers, are being used in the fore- 
casting of that component of sea-ice drift. The existing equations are 
modified and programmed for use with a computer to produce a forecast 
of wind drift, thus permitting rapidly made forecasts over large areas <> 
Applicability and limitations, ways of handling other factors, i.ec per- 
manent currents and melting and freezing, and methods for general im- 
provement of the programs are discussed. The programs developed are use- 
ful as research tools and as applied to forecasting. 
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1« Introduction. 



The problem of forecasting the wind drift of sea ice has been im- 
portant to man since he first attempted to navigate the frozen sea. Use 
of the polar sea for supply and commercial routes and as an operating 
area for nuclear submarines makes desirable faster methods of ice fore- 
casting for large areas. Zubov, Shuleikin and Fukutomi [1^ 2, 3], among 
otheiSjhave investigated the problem of wind drift of sea ice and have 
proposed various formulae for prediction of its movement o The U. S. Naval 
Oceanographic Office (formerly U. S. N. Hydrographic Office) has published 
ice-forecasting techniques, including those for forecasting wind drift, 
which are used by Navy meteorologists and others [4]. 

The Oceanographic Office method employs graphs and tables to obtain 
a forecast by hand calculation. Forecasting for any extensive area is a 
slow process and much subjectivity is involved. Data available for hand 
forecasts are rather sparse, making it difficult to do a complete forecast 
even if time were not a factor. 

A machine method of predicting the wind drift of sea ice would represent 
a significant improvement over present techniques; (1) it would rapidly 
provide a complete wind-drift forecast for an area; (2) continuity of fore- 
casts would result from retention of computed values which may be modified, 
in part, by new observations and carried forward into new forecasts; and 
(3) forecast wind information, to serve as input data, is available via 
numerical methods. 

This paper will present one approach to the numerical forecasting of 
wind drift of sea ice, together with an examination of its applicability 
and limitations. 
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2o Background. 

Many expeditions to the Arctic have contributed to available infor- 
mation on the wind drift of ice, geology, meteorology and oceanography 
of this area; these data have been used by later researchers. In September 
of 1893 Nansen's ship FRAM, built to endure the crushing pressures of ice 
floes, entered the pack-ice in latitude 78-50 N off the New Siberian Islands 
for the express purpose of gathering data on the motion of ice and the currents. 
In August of 1896 after drifting for three years with the ice and reaching 
85-55 N, farthest north for any surface ship, FRAM emerged from the ice 
into the Norwegian Sea east of Greenland completing its gathering of 
scientific data in the Arctic From 1937 to 1940 the Russian icebreaker 
SEDOV drifted over a route similar to that of FRAM„ In the intervening 
years several vessels had completed shorter drifts. 

During 1937 a four -man Russian team under Ivan Papanin drifted for 
nine months on an ice floe, the Papanin Floe; thus, the Russians pioneered 
in placing scientific stations on drifting ice islands [Dyson, 5]. Since 
1950 a number of these drifting stations have been occupied by the United 
States and by Russia. 

Zubov [1] utilized data obtained by the SEDOV expedition to obtain 
empirical rules for the wind drift of close ice (ice which covers most 
of the surface). His results compared remarkably well with those which 
Nansen derived from data obtained during the drift of FRAM. Zubov 
concluded that over a deep sea, away from the influence of land, the wind- 
caused ice drift is, on the average, at an angle of 29^ to the right of the 
wind and at a speed which is 1/50 of the speed of the wind. In his treatment 
of the drift of open ice (ice covering between two and six tenths of the 
surface) Zubov shows that the speed of drift is a function of concentration 

and increases with both decreasing concentration and increasing roughness. 
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Shuleikin [6], using dynamic relations and the drift of the Papanin 
Floe compared to winds^ derived a theory of ice drift in the open ocean. 

In a later writing [2] he examines his theory further by using data from 
the Papanin and NP-1 floes and obtains constants for his equations, 
Shuleikin's comparatively simple equations consider only ice fields of 10/ 10 
concentration. 

Since Zubov does not give equations for the speed of ice of varying 
concentration, we refer to Fukutomi [3] for a treatment of wind drift of 
sea ice under these conditions. Fukutomi *s formulae are considerably 
more complicated than those of Shuleikin, introducing stresses on the 
vertical surfaces of floes by wind and current and stress of wind on the 
exposed water; Shuleikin considered only stresses on the horizontal ice 
surfaces. This treatment by Fukutomi makes necessary auxiliary equations 
and more data on floe dimensions. Further, Fukutomi’s constants were 
determined from data obtained in the Sea of Okhotsk and thus differ from 
Shuleikin's constants, based on Arctic Ocean data. However, similarities 
exist in the theories and results. Fukutomi concludes that movement 
predicted by his equations agrees well with Nansen's observations. His 
results support Zubov in that drift speed increases with decreasing con- 
centration. Comparison of Fukutomi *s results for lO/lO concentration with 
those of Shuleikin shows rather good agreement, particularly for the angle 
of deviation of ice drift from wind direction Deviation angle is also 
shown to depend on concentration; this is a reasonable result, since an 
increased wind-induced current in exposed water areas directly affects 
ice movement. 

Review of the above papers and others indicates that the several 

theories and methods are mutually consistent; their different results are 

related largely to local differences in ice characteristics (e.g. roughness) 
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which are described by certain parameters* Therefore, the method developed 
here uses the equations of Shuleikin modified for concentration by reference 
to Fukutomio 
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3. Models for the wind drift and changes in concentration of sea ice. 

Shuleikin [2] has provided the basic equations required to obtain 
close ice movement. These are: 

tan = tan Q + B/u (1) , 



and u = v[^(§ cos a + cos g) 



k 6 
a a 

k 6 
w w 



( 2 ); 



where : 

\(r = angle between wind direction +180^ and the direction of 
ice drift 

u = ice drift speed 
V = wind speed 

k = coefficient of friction between air and water 
a 

k = coefficient of friction between ice and water 
w 

6 = density of air 
a 

6 = density of water 
w 

CD = angular velocity of the earth 
= latitude 
(jD = CD sin i> 

“2 

A = a dimensionless coefficient - 1.27 x 10 

N = n k^ 6 
a a 

^2 CD 6 

w 

§ = n k 

N CD ^/2” 

a = act n + 1) 

3 = n/4 - a 

m = mass of ice per unit surface area 

_ 2 

B = 2m CD (g cos g 4- cos g) 

k 6 cos B 
w w ^ 
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Figure 1 depicts the relationships among some of the above quantities. 

In deriving the above equation^ Sheleikin has assumed an ice concentration 
of lO/lO. He further assumed the ice drift is in equilibrium with a constant 
wind, i,e, an unaccelerated motion prevails, and all stress is applied to 
the horizontal surfaces of the ice, 

Fukutomi*s equations require a knowledge of: (1) total area of ice, 

(2) area of exposed sea, (3) plane surface area for each floe, (4) vertical 
area exposed to wind pressure for each floe, and (5) vertical area exposed 
to current for each floe. Moreover, one must be able to obtain or estimate 
the energy loss due to eddy currents. Most of these factors are not known 
by observation and are difficult to estimate with any degree of accuracy; 
therefore, it is not practicable to use Fukutomi’s equations idrectly. 

The equations are too long to reproduce here, and the reader is referred 
to Fukutomi [3] for further study. 

It is possible, however to approximate one set of Fukutomi *s curves 
for wind factor, u/V, versus wind speed, V, with Shuleikin’s equations by 
the introduction of a correction factor, E ~ |^1 + ( ) (1 4- 3V* x ,0lJ 

where c is ice concentration and V* is wind speed in knots, (c < lO/lO). 

The approximation is shown in figure 2 where the curves according to Fukutomi 
represent complete energy loss due to eddy currents; this corresponds well 
to that condition when area of ice is large compared to ice thickness. The 
resulting equations are : 

tan ^ = tan 3 + B/uE (3), 

and u = V[(§ cos a + cos 3) J E (4). 

w w ^ 

The modification was obtained in a purely arbitrary manner in order to best 
fit Fukutomi *s curves over the range of winds normally observed. 
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FIGURE 1’. Wind drift of sea ice force and velocity diagram(after 
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O WIND SPEED (m/sec) 

rIGURE 2: Comparison of wind factor vs wind speed according 

to Shuleikin and Knodle with similar curves according to 
Fukutomi . 



Now, a method is desired to obtain the rate of change of ice con- 
centration at a point and the resulting concentration at that point for 
some future time, Simpson [7] provides equations, which are readily 
adaptable to numerical methods, for computing these quantities. The 
equation of continuity of ice concentration given by Simpson is; 



be = 

at 



Su . bv . 




Sc 
""i Sy 



(5), 



and the equation for concentration at a given time is : 



where : 



be 

St 



= local change in concentration with time 



c = concentration 

u^= x-component of ice-drift velocity 

Vi" y-component of ice-drift velocity 

X = x-coordinate of location 
1 

y^= y-coordinate of location 

t.= initial time 
J 

tj_^j = time at which new concentration is computed 
— Sc 

t. = time over which — is computed, the time step. 

J St 



( 6 ); 



Equations (3), (4), (5) and (6) permit one to determine wind drift 
of sea ice and concentration changes resulting from that drift, thus 
completing our models. 



9 



4. Development of the computational programs. 



One of the two programs, ICEFORl, might be considered as a step to- 
ward the final program developed here, using the models developed in 
section 3. ICEFORl, reproduced in Appendix A, was written using equations 
(1) and (2) (Shuleikin*s equations) ; it produces wind-drift vectors 
for ice fields with k and k known and uses inputs of latitude and wind 
velocity. While it is not the final form of the wind-drift forecast, it 
does have some utility. 

The next step was to program equations (3) and (4), those modified 
for concentration, to obtain outputs of ice drift similar to the output 
from ICEFORl. Basically a program segment of ICEFORl was modified to 
accept all ice concentrations. Changes in concentration were next con- 
sidered. Equation (5) is put into finite difference form: 



— = time rate of change of concentration 
c = concentration 

u = X- component of drift velocity 
V = y-component of drift velocity 

D = distance between data points^^ and subscripts indicate data 
points as follows: 




(7), 



where : 



Ac 



2 . 



1 . 



.3 



0 



Solution of equation (7) requires values of ice-drift velocity 



and initial concentration, c , at each data point for each step. A new 
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concentration, c^ + Ac, can then be determined for each data point. One 

Ac 

notes, however, that in order to compute ^ at any one pointy, data from 
four other points are required. To satisfy this requirement, and to obtain 
an area forecast, a rectangular grid to approximate an area of the earth's 
surface was devised. The grid is a 546-point array of dimensions 21x26. 

It has an adjustable scale which permits the distance between grid points 
to be changed by use of an appropriate scale factor. Obviously if this 
plane is made too large it no longer approximates an area of the earth's 
spherical surface. However, for an area the size of Baffin Bay, requiring 
a scale factor of 2.0 for representing 30 nautical miles between grid 
points, the approximation to the earth's surface is good and no problems 
arise. The grid permits input of data over an area and further permits 
computation of convergence and advection within that area, as discussed 
in detail later. The program is set up to determine ice-drift vectors, 
changes in concentration and new concentration^ point *for-point throughout 
the grid; it is named ICDGRID, 

Several authors, notably Browne and Grary [9] and Fukutomi [3], 

have obtained values for the coefficients of friction (k and k in 

a w 

equations (2) and (4)) which differ from those used by Shuleikin. Also, 
Wittmann and Mac Dowell [4] apparently use different coefficients in modi- 
fying Shuleikin's equations for Baffin Bay. It is apparent that these 
coefficients do vary, as we might expect, for there are variations in the 
roughness of ice with age (newly-f ormed , thin ice is smoother than old 
ridged ice) and with location. Therefore, the programs developed here 
permit the user to enter values of these coefficients which best fit 
the area of forecast. 

Special problems arise when land boundaries are approached,* therefore 



11 



the program is written to recognize points of concentration entered as 
2.0 as land points and to ignore these points in computing ice drift. 

Certain of the land points are required j, however ^ in the finite differencing 
process used to obtain Sc/St o In order to treat concentrations near the 
boundaries the program is written to see the concentration of land points 
as lO/lO when the drift is onshore and as 0/ 10 when the drift has an 
offshore component. In this manner changes in concentration along land 
boundaries are made realistic. 

Program ICEGRID is outlined in Appendix B. It utilizes inputs of 
ice concentration, ice thickness and mean daily surface-wind velocity to 
obtain a five-day wind-drift ice forecast in one-day steps. The outputs 
are grids of ice concentration, direction of movement, and 24-hour distance 
of movement for each of the five days of the forecast. 
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5. Applicability 

The programs developed here are applicable to any geographical area* 

If properly indexed in latitude^ ICEGRIB will compute over the geographical 
area desired, and ICEFORl computes for geographical points as directed. 

Accurate forecasting in each case requixes that the coefficients of 
friction, k and k , which tend to vary with location, be known. Then, 
given primary inputs of ice concentration, ice thickness and wind velocity, 
proper outputs may be produced. 

There are t&^o basic uses of the programs: forecasting the wind-drift 

component of sea-ice drift and research. Both ICEGRID and ICEFORl provide 
forecasts of the wind drift of sea ice, with ICEGRID providing the most 
complete forecast. ICDGRID provides an area forecast and gives, in addition 
to direction and distance of movement, a concentration field: ICEFORl produces 
only the ice-drift steering vectors at selected geographical points. 

With regard to research ICEFORl may be used to obtain paired values 

of k and k for an area in which sufficient observational data of actual 
a w 

drift are available for comparison with computations. The rapid computa- 
tional ability available in machine methods makes this an extremely advantageous 
tool. Additionally, one may easily extend the possibilities for research 
with both programs to include investigation of permanent currents in 
the frozen sea, to study effects of specific features (i.e. pressure ridges, 
floe size) on the wind drift, and to relate drift to geographical features. 

The fact that ICEGRID is applicable to all concentrations permits its use 
in more general research than ICEFORl, which is restricted to concentra- 
tions of 10/ 10. 
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6. Limitations. 



This section discusses the limitations of ICEGRID as a forecasting 
method. Since only a limited portion of the total forecasting problem 
is treated in the theory used here, only that part of the drift directly 
due to the wind is found; components due to permanent currents and concen- 
tration changes due to melting and freezing are ignored. Some of the 
limitations discussed are, therefore, not related to the stated problem, 
that of forecasting the wind drift of sea ice. 

Permanent currents . 

Probably the most important limitation of the program is the failure 
to consider the effects of permanent currents. Wherever such currents 
exist they must affect the total drift of the ice. Therefore, the program 
should be expanded to include these effects. 

One may, of course, subjectively superimpose these effects upon the 
computer output, drawing upon a knowledge of the permanent currents to 
infer the total ice movement. A more direct method of considering these 
currents would be to input permanent current vectors at each grid point, 
adding them vectorially to the ice drift due to wind. The total ice 
drift and the resulting changes in concentration would then be given by 
the computer output. 

Curiously, permanent current data are not always readily available. 
Most current observations are made during periods when little or no ice 
is present, and these would likely be poor approximations for times of 
the year when considerable ice cover is present. Still, in many areas, 
say in the East Greenland Current or in the Labrador Current area, a 
calculation based on any approximation of the permanent currents would 
surely be of more use than one considering wind -drift alone. 
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Melting and freezing. 

A problem perhaps equal to that of permanent currents is that of 
assessing the influence of melting or freezing during the forecast period. 
Assuming no effect may be reasonable early in the melting season; however, 
it is most unrealistic at the height of the melting season and early in 
the freezing season. At these times, ice boundaries may actually move 
opposite to the drift if these effects are large enougho This time of 
year from greatest melting to early freezing is also the most important 
to the forecaster? as this is the time when ships ply the ice covered 
sea and the need for an accurate forecast is greatesto Although the 
problems are complex, models which account for ice melting and freezing 
exist and the necessary meteorological and oceanographic data may be 
extrapolated or approximated from climatology. A machine method of ac- 
cumulating degree days of frost or warming degree days might then lead to 
machine-produced estimates of freezing or melting to be incorporated into 
the ice forecast; thus a more complete forecast would result o 

Hummocking and ridging o 

The effects of additional hummocking and ridging are neglected once 
k and k are chose; actually, the drift characteristics of the ice are 
altered through changes in the roughness parameters. The immediate re- 
sult of failing to consider the above factors is a forecast of slower 
than actual drift in areas characterized by increased roughness; for, 
as previously stated, Zubov [1] concludes that increased hummocking increases 
drift speed o In Fukutomi*s development [3] hummocking and ridging are con- 
sidered in that they increase vertical areas exposed to both wind pressure 
and currents o Clearly, these factors should be given some consideration; 
but with presently available data any modification would be largely sub- 



jective, as aerial observations indicate only the presence and not the 
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dimensions of such features; further, one may not know the extent of these 

features present in the data from which the determination of k and k 
^ aw 

was made . 

It is also highly desirable to be able to predict the formation of 
ridges and hummocks, as ice dominated by these forms is difficult, and per- 
haps impossible, to navigate. Large hummocks and ridges also pose a hazard 
to submarine navigation by restricting water depth. They may lend to in- 
creased formation of fast ice, if they ground in shallow water, and thus 
alter the effective land boundary. This problem requires further study 
directed toward a dynamical solution, but will not be considered here. 

Transport from outside the area. 

No consideration is made of ice advected into the forecast area from 
beyond its boundaries, since no means is provided for entering data from 
outside the area into the machine. This inability to compute new con- 
centrations at the grid boundary points leads to corresponding errors in 
concentration near the boundaries. The problem may be further complicated 
by iterated use of the output as continuity data for new forecasts. If 
boundary points are not modified by new observational data anomalous fore- 
casts near the boundary can result, since the original input concentrations 
will have been retained at the boundary points. Accordingly, the numerical 
forecaster must be alert to change such points from climatology when such 
changes are required in the interest of continuity and common sense. Of 
course, one might use a grid scale to forecast over an area larger than 
that actually required; thus, the boundary points would have less effect 
on the desired forecast area. 

Nearshore drift. 

In the nearshore region of a water body, including that portion near 
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the edge of the fast ice, the drift of sea ice may be quite unlike that 
in the open ocean. Zubov [1] discusses the nearshore drift with respect 
to the effect of winds in piling up or driving away the ice from the shores. 
Shuleikin compares ice drift along East Greenland with that in the open 
Arctic and shows the differences which exist. From these and other dis- 
cussions, it is apparent that the presence of a solid boundary causes 
local changes in concentration and alters the direction and speed of ice 
driftc 

The effect of a solid boundary on concentration is to act as a nucleus 
for increasing concentration when there are onshore winds; and, contrarily, 
to act as a zero concentration point in the case of offshore winds (there 
is no source of sea ice on the land to contribute to offshore advection). 
This portion of the problem is considered in ICEGRID by forcing the computer 
to see concentrations at adjacent land points as 0/ 10 and lO/lO for offshore 
and onshore winds, respectively, in computing changes in concentration at 
nearshore grid points. 

The influence of the boundaries on the drift vector itself is not con- 
sideredo To make such consideration would .require that the program have 
available not only the location of land points but the orientation of the 
boundary as wellc To insert such information into the program would have 
required tying it to a specific geographical location which would have made 
it less general; therefore, no attempt was made to handle this effect. 

Thickness . 

The input of ice thickness is as much a problem as considering hummocks 
quantitatively. Few observations of thickness are available, from an oc- 
casional ship report during the navigable season or from a land station; 
these stations may report regularly on fast ice in the area. The result 
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is that thickness must be input from climatology or long-range ice outlooks 
formed from climatology and heat budget considerations » These provide 
reasonably good data over large areas ^ but local anomalies may be entirely 
smoothed out, especially by the use of climatology « Fortunately my compu- 
tations using Shuleikin's equations show that drift distance decreases 
on the order of 0.5 nautical miles per day as ice thickness increases from 5 
to 8 feet, for all wind speeds o The angle of deviation from the wind increas- 
es with increasing thickness; however, for a thickness change from 5 to 8 feet 
the maximum difference in the angle is about 15^ at low wind speeds and less 
with greater wind speeds, being on the order of 5^ for 20-knot winds » Thus, 
the climatology or ice outlook thicknesses will give generally acceptable 
results when used in computations o 

Transient states in the ice fieldo 

In the equations for sea ice drift it is assumed that the ice moves 
everywhere in equilibrium with the mean wind over the period of a forecast 
step, in this instance one day. It should be apparent that, due to inertia 
of the ice floes, a transient stage occurs whenever a change in wind velocity 
occurs » Zubov, however, states that, ”all observers are struck by the speed 
with which ice reacts to any change in the wind,”^ indicating transient 
states in the ice movement to be rather short. They are ignored, therefore, 
since the program steps are of one=day duration. It might be necessary to 
consider these states for shorter program steps. 

Observational data. 

Any forecasting program may be severely limited by the availability, 
accuracy and form of data to be used as input. Currently available obser- 
vations of sea ice limit the forecast in all of these respects. The majority 
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Zubov, N. N. , Arctic Ice. 



Moscow, 1953. p. 355. 
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of data are from aerial reconaissance flights made at somewhat irregular 
intervals over the limits of the ice; thus the amount of data available 
is adequate over only a small areao 

In addition to data scarcity, data reliability is often open to 
question, since the location of ice observed from the air is directly 
dependent on the accuracy with which the position of the aircraft can be 
fixed; and concentrations depend on the accuracy of the human eye. The 
Leningrad Division of the State Oceanographic Institute has tested an 
instrument designed to measure ice coverage. In the tests, during which 
aircraft position could be accurately determined, it was noted that in- 
dependent visual observers not only differed with observations made with 
the instrument, but among themselves as well. Errors in positioning the 
ice boundary were up to 15 miles; errors in area covered were up to 40^ 

[Betin; et al, 11]. The implication is that even direct observations 
may reflect the actual condition inaccurately. 

Current ice reporting messages are not readily adaptable to digital 
methods. A report must be plotted on a chart before data can be extracted 
for program inputs. This is a time-consuming process and results in a 
great loss in the speed advantage realized by a machine method. 

Programming language . 

The program as written is in the programming language called FORTRAN *60. 
This language is simple and easy to use and has speeded the production of 
the program. However, it has some disadvantage in its use; for example, 
all data must be input via punched card to magnetic tape, and computed 
values cannot readily be retained for use in new forecasts » Also, since 
a wind- forecasting program is not available in the language, all wind 
forecasts must be extracted from prognostic charts prior to being input. 
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All this results in loss of time in entering data into the program; the 
entire process requires more than 15 hours for a 600 by 750 nautical-mile 
grid, provided the required ice and weather charts are already plotted. 
Program step size. 

Program step size is only a limitation in that it affects other parts 

of the calculation. We lave mentioned briefly transient states of the ice 

field; it is possible that the wind field might be such as to make the 

transient state important even in the 24-hour steps now taken, and these 

stages surely would be important if the steps were shortened to the point 

that steady state did not apply. The step size is important in the com- 

^c 

putation of As now programmed this quantity is obtained by a for- 

ward differencing method which can introduce errors over a large number 
of time steps. If the number of time steps were increased, a significant 

error probably would result, necessitating a revision of the method of 

3c 

obtaining . A central-point differencing scheme is suggested as a 
means of accomplishing this if program steps are increased in number. 
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7. Results. 



This section discusses only ICEGRID as ICEFORl simply solves Shuleikin's 
equations. Due to the time involved in extracting data for the program input 
only three runs were made for verification. This small number of runs does 
not permit statistical analysis of the forecast versus the observed ice 
field. Analysis is further complicated by the fact that observations often 
cover only a small portion of the area of forecast; hence forecasts were 
made for the few periods when verification data were available. Despite 
these disadvantages, several items of significance are shown by the results. 
The effects of permanent currents, for example, are clearly demonstrated. 

First to be considered, however, is the degree to which we can depend on 
the program as a forecast of the wind-drift component of the sea ice motion. 

In the analysis of results it must be remembered that only the wind- 
induced component of drift is predicted. The immediate conclusion under 
this condition is that ordinarily one would not expect the forecast to 
conform exactly with the observed ice conditions. Such a conformance 
would result only if the drift due to permanent currents and heat budget 
effects were negligible. 

Using coefficients of friction obtained from ICEFORl, ICEGRID produces 

almost identically the drift which one obtains by hand using the curves 

of the U. S. N. Oceanographic Office [4] for Baffin Bay (these were for 

k = .003 and k = .013). Thus, one is assured that for areas of lO/lO 
a w 

concentration forecasts of wind drift will be at least as reliable as 
those obtained by hand. Considering that we alter Shuleikin's equations 
for concentration, we may extend this reasoning to say that forecasts 
of wind drift ought to be better than those obtained by hand since hand 
forecasts are made using curves for lO/lO concentration, without modification. 
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While no hand forecast has been made, for such comparison one machine fore- 
cast was done considering all concentrations to be lO/lO for drift calcu- 
lations; a significant drift difference was noted from a forecast modified 
for concentration, using the same input data. The forecast considering 
all concentrations was better for all points in the grid than the one 
assuming ice concentration to be 10/ 10 everywhere o 

To further establish the worth of ICEGRID in predicting wind drift, 
an examination of the machine forecast versus a forecast based on persis- 
tence of the input field was made. Figures 3, 4 and 5 show the forecast 
boundary versus input and observed boundaries for each of the three fore- 
casts made. In each forecast the machine result was better than persistence, 
in that it gave results which were generally more like the observed ice 
field so far as ice boundaries were concerned. The movement was generally 
in the right direction. Comparing forecast concentrations to persistence 
of concentration also shows the forecast to be better since it gives, at 
least, a hint as to the direction of change in concentration. Figures 6 
and 7 contrast an input concentration to its five-day forecast. 

In comparing the machine forecast with the observations at verifying 
time one finds, in general, a considerable discrepancy in the location of 
ice boundaries, say five to 30 miles per five days, and sometimes in con- 
centrations, as much as six tenths, which must arise out of effects other 
than pure wind drift. In most instances, the forecast drift is of the right 
sense, but it is generally insufficient to produce the observed results. 

A full analysis of the forecast results versus observation requires, then, 
a consideration of factors ignored in the development of the method Most 
of these have been previously discussed (see section 6) and will, therefore, 
be mentioned here only as they seem to affect the validity of a forecast 



made using only wind drift. 
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The effect of permanent currents is the most apparent of all effects 
exterior to the wind. Forecasts generally show the observed ice boundary 
to vary from that forecast by an increment which, for Baffin Bay, appears 
to conform well with the permanent current believed to exist. It is con- 
ceivable that for extremely light winds the effects of the permenant 
current would dominate leading to a drift entirely different from that 
forecast o In running forecasts, it was noticed that in central Baffin 
Bay around 67 N latitude forecasts conformed rather well with observations, 
indicating a lack of permanent current there, or at least a current parallel 
to the ice boundary. The study of permanent currents through use of this 
program is suggested as a possible research problem. 

A second noticable effect is that due to melting, since the forecasts 
were all made during the melting season. This is seen most easily through 
comparing observed concentrations with those forecast. Large decreases 
beyond those forecast are sometimes observed Such changes are likely 
due to melting and not to divergence or advection alone, as no reasonable 
amount of drift could account for them. 

In a previous section data errors were discussed. These contribute 
to apparent errors in the forecast. Since there is, at present, no method 
available to exclude such data errors, one must use all data assuming it 
to be correct. Fortunately, with regard to ice concentration, it is the 
navigability of an area which is important and not the absolute ice con- 
centration. Thus, some error in forecast concentration may be tolerated. 

The forecast concentration has value if it produces an indication as to 
the navigability of the frozen sea. Analysis of program outputs indicate 
the method to be capable of such indication. In the forecasts produced 
it has always been possible to determine rather accurately areas of de- 
creasing concentre tion* hence, one is able to determine the areas most likely 



to be navigable. 
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Analysis of forecasts indicates the input wind fields to be poorly 
defined with respect to divergence and curl of the wind on the scales 
desired. Thus, small-scale features in the ice field are not shown or are 
smoothed out over the period of forecast. These small features are often 
transient and highly important in the local navigation problemo This raises 
a question as to the effect of small-scale features in the wind field 
on the movement of sea ice. Browne and Crary [9] indicate that the ice 
fields may move with the average wind over a 60-nautical-mile square 
surrounding the location, but are not able to determine conclusively 
this relation, as weather maps fail to supply adequate wind information 
on this scale; and, moreover, all sea ice is not in the form of ice 
fields. It may be, then, that wind forecasts to include sma Her- sea le 
features in the wind field will be necessary before the smaller- scale 
features of ice drift will become apparent. Since weather reports in 
the Arctic areas are few in number such a definition of the wind field 
may never be possible. 

The results obtained make ^-t apparent that areas of Baffin Bay exist 
wherein nearly all the drift of ice is that due to wind. This seems par- 
ticularly true in the central area near and just north of the Davis Strait. 

In such an area the wind-drift forecast should constitute a total forecast 
if no heat budget effects are present. 

The few forecasts made using program ICEGRID indicate that the method 
produces good results for wind drift and that it is possible to expand it 
to include other effects. 
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8. Recommendations . 



A review of the preceding sections invites consideration of ways of 
improving the program. One is struck by the need for better methods of 
obtaining and reporting datao Additionally, results obtained point out 
the necessity of expanding the method to provide a complete ice forecast. 

A study of methods of observation and reporting of data has not been 
attempted here and might be the subject of separate investigation. A 
few possible aids to improvement of data quality and availability are 
worthy of note, however. Stereographic aerial photography and TIROS 
(Weather Satallite) photographs are suggested as providing more complete 
observational data„ Stereo photographs would give detail as to the dimensions 
of features in the ice such as hummocks, pressure ridges, and even an estimate 
of thickness where ice and open water lie adjacent. Large scale concentra- 
tion data should be available from TIROS photographs where cloud cover does 
not obscure the surface. Scanning of photographs might be done electronically, 
obtaining directly concentration inputs for each grid point. 

Even with better methods of obtaining data, time required to input 

that data to ICEGRID will remain large if the reporting format is not 

changed. The present format presents location, amounts and characteristics 

of the ice verbally, requiring the data to be plotted, extracted for each 

grid point and punched onto data cards prior to computer input. A suitable 

format, from which a direct input to the computer could be obtained, might 

be similar to that of the World Meteorological Organization (WMO) used to 

report weather. Such a report format is suggested in figure 8 and could 

be included in the WMO codec This format would report total concentration, 

concentrations by size and boundary location points for the concentration, 

all reported digitally. Additional plain language data could be included 

at the end of such a message. The code groups could be repeated for each 
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C Cn, P2n3 ULaUUt L.LaLoUt 

C C .* concentration in tenths 

! tenths of slush, brash, block 
rij tenths of small to medium floes 
p 3 ; tenths of large floes field 
UU ! latitude in whole degrees 

UU: longitude in whole degrees as 
in World Meteorology code 

t ‘tenths of degrees indicator as 
below: 
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FIGURE 8: Suggested Ice reporting format 
to be used with the World 
. Meteorology code. 
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"concentration area" of ice observed^ The boundary between ice and open 
water could be described by similar latitude- longitude groups following 
a code group, LIMIT. 

We must also seek to provide a means for machine retention of com- 
puted data and a wind -forecasting program. Fortunately, the U. S. Navy 
Fleet Numerical Weather Facility has available a program computing a wind 
prognosiSc This suggests rewriting ICEGRID in SCRAP, the programming 
language used by FNWF, obtaining immediately a machine forecast of winds 
and, at the same time, the ability to retain computed data as continuity 
for new forecasts. 

Extension of the program to include the effects of permanent currents, 
heat budget effects, land and fast-ice effects on the drift vector (in 
short, a complete forecast)is also suggested. Such expansion would require 
considerable time and effort on programming and on obtaining the necessary 
topography and oceanographic datao The production of an efficient and 
useful forecast method dictates the need for this effort, however, and en- 
courages that time be spent on such expansions. 
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9. Conclusions. 



In the foregoing sections the programs presented here have been 
examined and some of the factors foreign to pure wind drift of sea ice, 
which one ought to consider in a complete forecast, have been discussed. 

The programs developed here are concluded to be superior to hand fore- 
casts on speed and on accuracy, due to the modification made for concentra- 
tion and the ability to retain computed values to provide continuity 
to the forecasts. With the inclusion of the effects of permanent currents 
and heat budget, one may expect to obtain a complete and useful ice-drift 
forecast, provided good inputs are available. 

Further, it is concluded that the ability to obtain values rapidly for 
the coefficients of friction, and k^, leads to the programs* great value 
in research work. 
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APPENDIX A 



PROGRAM ICEFORl 

Description . 

ICEFORl utilizes Shuleikin*s equations, without modification, to 
compute the wind drift of sea ice. Provision is made, however, for entry 
of coefficients of friction suitable to the area of computation^ Drift 
vectors are computed for one or more locations provided the primary in- 
puts of latitude, ice thickness, wind speed and wind direction are 
entered. Secondary inputs of longitude and date/time group may also be 
entered for reference purposes, but these are not required in any cal- 
culationo 

This program has its primary use as a research tool. Given a set of 
observations of ice drift, corrected for permanent currents, and of corres- 
ponding wind velocities, one may find the Shuleikin curve with the pair 
of coefficients of friction, k^ and k^, which best fit the data. Having 
obtained the constants, one may obtain ice-drift vectors to be used in 
forecasting . 

Table A-I defines the variables used in the program language. Where 
variables are equivalent to those in Shuleikin*s equations S)nnbols are 
included. The complete program is listed following this discussion. 

Data input . 

Prior to use of the program, it is necessary to assign values to the 
constants AKA and WKW in the subprogram, DRIFOR. The values assigned these 
constants in Figure A-1 are for Baffin Bay. The values assigned to RHOAIR 
and DSUBW may also be changed, but this is not usually required. 

Data input is accomplished using IBM cards with punched values of 
variables entered on them. The first card specifies latitude, longitude, 
date/time group, ice thickness and wind velocity, assumed constant for three 
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hours c The second and subsequent cards are available for specification 
of additional three-hourly wind data at the previous data point. Follow- 
ing the last card specifying wind data is a card on which wind speed is 
entered as 5000. This signals the computer to proceed to the next card 
specifying the complete data, ioeot another ’’first card”. When all data 
have been read, a card which specified latitude as 300 signals the computer 
to read no more data cards o 
Output of results 

The program prints output in columns, titled as follows: 



Result 



Units 



latitude 
longitude 
date -time 
ice speed 



degrees 

degrees 

day-hour 

knots 



ice direction 
distance/3 hrs ^ 
psi 



degrees 

nautical miles/day 
radians 
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TABLE A-I 



Definitions of variables in ICEFORl 



PLAT 

PLONG 

DTG 

THICK 

VELWIND 

DIR WIND 

VELICE 

DIRICE 

DIST3 

AKA 

WKW 

RHOAIR 

DSUBW 

ENN 

C^IEGAB 

XI 

ALFA 

BETA 

SQMASS 

BEE 

PSIO 

PS I 

WICE 



= latitude 
= longitude 
= data and time group 
= ice thickness 
= wind speed in knots 
= wind direction 
= ice-drift speed in knots 
= ice-drift direction 

= 3-hour drift distance in nautical mile 
= coefficient of friction between air and ice, k 
= coefficient of friction between ice and water, k 
= air density, 6^ 

= water density, 6 
= N _ ^ 

= O) sin = CD 

= § 

= a 

= B 

= M 
= B 

= an angle, used in computation 
= angle between wind and ice-drift direction, t 
- ice-drift speed in meters per second, u 
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APPENDIX B 



PROGRAM ICEGRID 

Description , 

ICEGRID computes the wind drift of sea ice using equations after 
Shuleikin, as modified for concentration less than lO/lO, and concentra- 
tion changes according to the equations of Simpson. Drift vectors are 
computed over a 546^ point array using latitude, ice thickness, ice con- 
centration, wind speed and wind direction as inputs. The purpose of the 
program is to obtain a five-day forecast of the wind drift of sea ice. 

Table B-I defines the variables used in the program language. Where 
variables are equivalent to those of Shuleikin or Simpson, symbols are 
included. The complete program is listed following this discussion. 

Data input 

Data input is accomplished using IBM cards with the values of the 
variables punched on them. Four constants must also be entered into the 
program; these are also input via IBM card. 

The first data card specifies the constants, SCALE, which is the 
scale of the grid, and BLAT, the base latitude. A scale of one provides 
grid squares of 60 nautical miles on a side; a scale of two gives grid 
squares 30 nautical miles on a side; and so on. BLAT establishes the 
northern-most boundary of the grid and is always the latitude of that 
boundary plus the inverse of the scale factor. The second data card 
specifies the values assigned to the coefficients of friction, AKA and WKW. 
The fourteen data cards following the two just defined specify the concen- 
tration of ice at all of the 546 grid points, values determined by 
averaging observations over the four adjacent grid squares. Land points 
are entered as 2.0. Following these cards are fourteen cards specifying 

ice thickness over the grid, taken to be the mean climatological thickness 
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for the period of forecasts Thickness is entered in feet and is obtained 
in the same manner as concentration for each grid point. Fourteen more 
cards specify wind speed in knots, obtained from prognostic charts, for 
each grid point by averaging over the four adjacent grid squares. Fol- 
lowing this group of cards is a group of 21 cards specifying the direction 
toward which the wind is blowing. 

The remainder of the input cards consists of groups of 14 and 21 cards 
respectively for wind speed and direction for each of the remaining fore- 
cast days. 

Output of results 

Results are printed in grid form with the outer boundary points of 
the grid omitted. The input concentration grid is reproduced as the first 
output. This is followed by grids of concentration, direction of movement 
and distance for each of the five forecast days. 
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TABLE B-I 



Definitions of variables in ICEGRID 

CONCICE 
VELWIND 
DIRWIND 
THICK 
TEMPICE 
VELX 
VELX 
DIST 24 
SCALE 
BLAT 
PLAT 
VELICE 
DIRICE 
DCDT 

YICE, YYICE, SICE, XXICE = temporary concentration storages 
AKA - coefficient of friction between air and water, k 

a 

WKW = coefficient of friction between ice and water, k 
RHOAIR = air density, 63 
DSUBW water density, 6 
OMEGAB = CO sin 0 , CD ^ 

XI = § 

ENN = N 

ALFA = a 

BETA = 0 

SQMASS = m 

BEE = B 

PSIO = an angle , used for computation 

PSI = angle between wind and ice direction, ^ 

CONFAC = factor modifying for concentration 
VVICE = ice drift speed in meters per second, u 

C = constant used to modify outputs of t|f and u, if consistently 

low or high, normally set to one 



- ice concentration, c 

= wind speed in knots^ V* 

= wind direction +180 , in degrees 
= ice thickness, enters into computation of mass 
= a temporary storage for new concentration 

- x-component of ice speedy u, in knots 

- y-component of ice speed, r, in knots 

= distance of ice drift in 24 hours, in nautical miles 
= scale factor, S, determines grid size 

- base latitude permits proper orientation of grid in latitude 
= latitude , deg ree^ 0 

“ ice drift speed in knots 
= ice drift direction, in degrees 

- change in concentration with time, Sc/St 



42 



II H jf ?f 




U.H- 

U=) 

0 . 

a«-o 

M *-«LU 
K tl-X 

Q^OOUJ 
►-<oc: • 
zo^x 
Luujax 



oc 

u 



00 

< 

o 



o*^ 

i-<CNJ 
X •• 
h-O 



OQi 


UJX 








ZOOQ 


oco 


— J- 






OU.ZZ 


Oh- 


CNJ CNJ 






O << 


u» 


*>- 






>- 


z 


OtO 






Z<OH 


UJm 


<M*-I 








UN- 


—o 






UJO 




o 






mmci-o) 


X 


z— 






ui to 


o 


. 






o< :3 


oc 


:xcNj 






Z OiX • 


u. 


QC • 






<QCZXUJ 




I -.0 






XO*-#H -0 


X 


OCM 






OU.X < 


H* 


•iw* 






^oo 


X 








omtoujto 


o 


— *-J 






zamo ^ 


to 


CNJLU 




o 


<UJUJ <0 








CNJ 


H-ZOZ 


o 


o •. 




•k 




z 


CNJ — 






oc o to 


< 






II 


OC 5 C«-Ht 03 




OCNJ 






H-XX'-i 


K 


z •• 




•k 


O H- O 


to 








ujjs’ vutum 


< 


XCNI 






>CNJOQc^m 


UJ 


-J^ 




CNJ 






UJX 




9» 


oz<kk-n. 


o 


>-J 




r— • 


Z*M - 




••UJ 




II z 


zi-<o • 


QC 


—3 




•JO 


Qcmo< «to 


Oo 






•k— 


UJH*^OUJ*Z 


r— 


CM— 


• 


— h- 


UJZH 0CO< 


XUJ 


> 


o 


o< 


1 — •<z»-» 


oo 


Z OCNJ 




»o: 


tOOQCtO <0 


CMnh 


O CNJ •» 


no: 


•-iH- 


0 .h-> 00 < 


•-X 




<o 


wZ 



O 

<NJ 



^00 

(NJOO 

•iUJ 



to 

CNJ 

•i 

CM 

M 

£ 



O 

CNJ 

•» 

CNJ 

11 

z 



Z-- 



o 

CNJ 



CNJ 



II 



o 

CNJ 



CNJ 



II 



rO<^ 

mm 

•» «» 

CNJ.— 

mm 

•> •» 

CNjJt 

mm 



m ^<r^-iLiJry' 

UJOUJQCtO 

O II O X< 

*-»r 3 o:z»- 

c 3 ^acj>ocoujo xi-i 

ox UUXZOO V.K 
UJOO U.H-<tUJ 
UO^—UJ OIU V*< 
^ LLOmULUJCr vz 
X«-«««» OOO^N* 
X<0 U. wUJfOw 
<OCtU^CJO 0 UJCJ h- 
o:oQ.< z-j 
OOoO»-»XOOOOZX 
OQl l-0«-«ZH*^C£> 
O^OLH-t— < 00 <C^ 0 CO<D 
a <ZLUZtfUJCLU. 
m *-» iij > 

•-•UJ OCX z 
XUU.O(-* O 
H-MOU.O u 

uuouu u 



< 

:za: 

Ol- 

CVJZ 

»-UJ 

V.Z 



OLICNJ ^ 
u«-' »OJ 



X’v riC oo 

» • 9 

UJX. -> -) r-CM 

UCNJ — wZ I I O 

tuuj t-ity i-i • Q oo • 

X oo *^0 O'O 2 Zi-I i-ii^f— 

*-<UJ ^X ^ *-<Z iiC< 2UL »-• • H 

oo uj»- ^3c-*oa^oo::--‘aoomxo— XO-* 

m #fcmzo^*-<ujoo^ •^-jtuoocujoo*— 

OQ-^<0C •< »0 »X> »w a^UJUJ •»«<ClC tCNJCMTOS^ 

oxcNjoo o:^ooocNjh-<cNjwCNiii >o.CNjQ«~im •» •ww. » 

UJ •-00LL •»— <iJLw LL •* »0 wOOU-w QU,f— r-LUUJi— I 

Zh-o • uju. cr— on moU*^ 

O'VO— *0 ^CM— oiocNJCNjCNj »h-ro •*<rr)r^LL ••oro «ia(Nji-4i^Mi~.^ 
Ow Nfc»-i^wOOZ5^0^0>-<w»— CJwfww^CNjZ*^C02^*0ic:OOUJ 
H* h* ►— ZZO 

►-< ^ZCNJO^O*-'<.— < fc-H< j/^<K-< 0 ' x< X<^^ 00 »-| 

zxh-s^uj c3<HXoxo xo xzx*-o xo xmmooo 
*^oC 3 xooc<LL<oc<QC<za:<za:«-'c^ <zc^< 2 oc: wwz 

OCOCL i-^CNI'-iUJ -JOUJOUJ*-*OiiJ»-'OC^OOliJ*-^OUJ»-»OOCDlL.LLa 
aiLz Owooruj u.c^u.c^ u.o: u.atLOtt: uloo^h^o 

»-CM f-CNJ 



fO 



o 

o 

o 



OUJ 

<to— 


CNJ 


o 

<o 


o 


r- 


o o 

<CMO<fO 


o<o 


o 


LU^ 


Ujr- 




UJr— CX)UJ— 


to coo 


o 


CC 


QC 


N- 


CC CNIQC 


uu 




u 


u 




u o 



ro^ 

mm 



43 



•> 





• 




2 












t-. 




> 




< 
















u 




a 






• 










lU 




2 
















> 




» LU 




















-- a. 
















H- 




2 LU 
















u 




2— O 
















UJ 




•*3: 
















Q. 




X2 in 




















X JC K 












* 




LU 




w ► 2 
















QC 




•- 

ox o 




















»-•< .a, 










j 






•k 




X •• 










1 






o 




h-— o 




















•*2 2 
















o 




^2 < 
















H- 




2 •• J 

2X 
















H- 




••X >- 
















LU 




Xw < 
















m 




XJ* 
















UJ 




wCM Q 
















oc 




Oh- . .LU 

2 in in 
















•» 




003 
















r“ 




22 

ct: • i-#— luj 






' 










2 




•-I-. X3CQ 














1 


< 




02»-Q:a: 
















X 




•*2 UL ^ O 
















h- 




^ •••-•QQI- 




•0 


0 




in 


0 








2Xo:ww 




h- 


00 




CO 






q: 




2XQUJUJ2 




•» 


•k 




•k 


•> 




UJ 




-J-iO 




m 


f— 




0 


r— 




►- 




XLULUOO*— ' 


0 


fw 


00 




CO 


c^ 




< 




X0022I- 


»► 




•» 




•k 


•» 


r— 


LU 






f^ 


0 


0 




in 


0 


in 


a: 




oc3i: 0 : 0 : 0 ^ 


•» 




00 




00 


C7> 




o 




2‘^LU*—’—l— 


0 












r— 


•» 




•-iO^-U-U-2 




000 


OrO 




oco 


om 


in 


o 




:x •►r)2inLU • 




•r- 


• CO 




•00 


• C7^ 








_jajQ.i«iOOLn 


0 


CM 


—CM •* 


-J 


CM •» 


—CM •. 


o 


2 




ujoxinoxo 


• 


1 


•J 1 ro 




1 00 


•- 1 CM 


in 


< 


UJ 


*02 


CM 




^—00 




— *00 


1 — C7^ 


o 


X 


u 


•‘-JOLULUU*-* 


1 


•j •» 






f" «k 


-J»— •* 


0^ • 
•—o 


l— 


<r 

o 


^>o*-^'^in 


G 


Xf:: 


^ VCM 
1 *00 




* fe 
-jeo 


v' 4-CM 
•-JO' 


M 


m 


m 








-)f— — 




id:— 


->x— 




(no»- 


V. 


a-- Luuj2a:oo •» 


1 


1-4 ^ 


w 


•k— <k 


«k4^ 




Lucgcg 


XUJ^^X>>*-OCMCVJ-> 




w-n J 


LU 


-i-j 






•J 


XOCM2< II H XX •* 




♦-•id: 




u 




LU^-i^ 






x^q: - 0 ^-— a: in cm eg w 


•» 


Uw » 




w »• 


0 — • 


UJi-« 


to M II 


1 i-as:ozzuju. n 


II LU 


LU-> 


i-iUJ-J 


u 


LU-> 






231.2 


h- »-i U. X u 2 2 1- IX -J O 




OU«-*CJ 


CJ 2 


o>-t 


00— 


l-Hwl 


JJOX2 


<h-UJwC3^ » •»LUa*-:^*-'OCM»-'wO 


0 2^-*^ ♦ 


• 0 


l-IwO 


0 2^^-^ 


OUJ^I-I 


J<OLUa.XXQ^ 


0 


• 00> •C7' 


•c^ou>•- 


aa^u 


ox •CT 


•O'OOX 



icoiCK-'-^aetD-ji-'oaaa.a; ujcM<Mzt>^2_j^r-or-oz-< ii cu ii cu ii , z-Jt>03»-aDOZ_J 
0>— ►-*<005: II I— rj— 'WOQCOOO H OUJH II M OUJUJ UJ UJ'OUJII M II olu 

OOI-Ot^»- II KU5_JO(/>X>-ZO--'-0l-OU>lUOUJ0UJU>0O0OO 00UJ0UJ0UJO3 
*-ZZ|- Z<UJ_JZ wO|-wwOI-UI-O—-*-'l-'-il->-i!ww0l-OI-Oww 

U.OOZOOZ_ll-<O(/)UJUJ3;00OOU.OOLLU.i->OMa>->u.U.>O>O>-|U.U.MOMO<->U.U. 
•-•OOllJOOX0LOUOa:3>OZOO'-iQO>-»<-i>C»-C»'«"-i>O>C3>- •-“-'XOXOX— <-i 



0 


X 


•-UJ 


-JO 










i 




CMO «-2 


X 


>- 


•— -J 


h- ^inoo 


r- 


* 00 '— CM 




0 Jt*OCO 


inc'f- 


minmo 


0 


2 


002 




h- 


f-h-CDCO 


00 


QD, 000000 


00 ooooc^ 


0 


0 


UJ 


wllo 










1 




} 0 


u 


u 


00 















44 



o 

UJ 

QC 



UJ 

O 

• K 

OX 
<NJ I 
r-LU 

-^x 

4-«-- 



UJ 

U 



UJ 

> 

UJ 

o 



a 



a 

oi 



a 

o 

z: 

I/) 

z 

o 



< 

o 

o 

UJ 



*■ 


















0 


CO 






co 
















► 


z 




>1 


-JX * z 
















X 


UJ 






1-4 •‘UJ 0 
















u 


X 






— »-4 
















►•4 


»i-4 




\ 


Xl-4< h- 
















X 


UJ 






-1—0 < 
















1- 


-1 




j 


UJXOO QC 
















•» 


0 




• \ 


0-J« H 






— 










Q 


X 




i 


1 UJ-- KZ 






in 




— 






z 


CO 




i 


OUJ 






CM 




in 




in 


I— • 






3 


1 • 00 






► 




CM 




CM 


3C 


0 




1 


1 ---0 oz 






CM 




► 




• 


oc 


z 




1 


-J «CM *0 






M 




CM 




CM 




4-4 




4 

i 


«o 






»-i4 




II 




II 


0 


00 




■a 


XM’N. Jd- 






► 




X 




Z 


•k 


0 




.1 


CM3C 






*— 




► 




► 


0 




• 


] 


•-'VUJ -f UJ 






0 




— » 




— * 


Z 


UJ 


0 


I ' 


—0 •-«z. 


0^ 








0 




0 


1-4 


u. 


zm 


1 


X— Z-J 


z 




•• 




CM 




CM 


3 : 


1-4 


Oin 


\ 


-JJ>0^^ 


►i4 




CM 




•» 




•k 


-1 


a: 




1 • 


UJ^>*-4 ►UJ 


► 




II 




CM 




CM 


UJ 


0 


r-ujh- 


1 


. 0 ► 1 i--)o: 


z 








II 




II 


> 




h-OO*— 


1 


wf-UI<l-40 


1-4 




► 




•J 




Z 


*> 


UJ 


^ 0 • 




— 40q:— K 






— 


0 


•» 






1- 


0 


•!/) 


j 


— « -7»-4H-UJl/) 


UJ 






-J 










—•.4 


inuj 4« z 




* 


0 




► 


UJ 


-1 




Z'vv. 


-13 


•a.o< 






f-H • 




•-H 










Q-XCO 


» z*-* 




-J>^>— UO-J 


aoo 


0 


w 


u. 


X 








301 


o:co<to 




X *--1 ♦ ZZ*“» 


XI- 


•J 


UJ 




W 


\ 


— CM^ 


CM 


►H- 


Xq:o< 




►-nuj-*ono/ 


ujo 


UJ 


1 ) 




UJ 


0 


• • 


o:<r*n 


5TUJUJO/ 




•^»-4>-jOu<i 


HCL 


►-« 


»->4 


h- 


0 


-J 


CMO'O 


Qxa. 


zf=o 




»-»— IX mo: 


II H- 


u. 


0 


Z 


H4 


Ul 


l-U-U. 


u.<x 


>UJ II o 




wuj— ►z— otno^o 




z 


r—UJ 


0 : 


r-* 


in 00 CO 


u 


►0 


• M S — K 



O xz: H-l XU. ►-•r-r- 

u ouj o o o 

w <Nl> W csl wCM*— 

w r— O ^ r—UJ W • • 

►CT'CT' 

a^^ox. >ou.u. 



uJO-J“^*-»-ja.CNj<Njzo X 

oz 

•-IO ••^UJ •.UJCJCM (NJO 

0Z^>-0->h- Mil XZ 

O O ZQ I ^Z*— XZ*-^‘-<C\J'vh- 

• • oo'^uj<wa»-»*-»wQccsjN^< 

0 ' 0 »->ou I *-<>xujz ujCLcgwa-^o^m— w<-.oww 

II ^ II f— H w.- I oo<irnrno »— i-u-*v h-i— v. i— h- 
UJ UJ UJ||>— OO*-*O^-<ZN.I-K<O'Vl-h-Ct0's*h-<< 

OOOOOH--J •</>aoo»~^ozzxuj^zzx v^zzsih^vzxx 

*.^v-i-ii-i-40ujouJXUJ Z*-'*--iQ!:o i-i^csio »-»a:a:QcDa:0ZKt/)Z»-0^00 

xoxoxu>c\jh-uji-ooo3ca:oz a:a:oz 0 : 0:0 o:aozo«-»i/)oxuj<o:o 11 Xoo 

xoxoxow^ 3 i. 30 ouoa.u-o O.O.U.O Q.au. GLU.u.uiu^o u>e)o:ujo:uo:o 



O^OZ-* O CNJ 
«-»<NJ 00 : 0 X 1 /)-- 

|-UJ‘-iXH<UJU. 
UJCOZh-X OUJ|~ 
2 —I CO UJ ct CX 
i-iai-o:>h-QOo 
K ►Ol-wOwUJUl 
OUJUZU.ZUJLJ II 
ouo:uj(/)b<:«j CM 



m 

I 

UJ 

mm 

• o 

II •- 
OC II 

•—x 



o:«-^cDoo: ocoi- •<cu 







•-cMa ql. 


-1 


•-CM 


—CM 




•CM 


f— CO 


z 


> 


CM 


m 


0^0 XCMX 


m-j 


CM 


m 




^a>o 




< 


z 




O' 


c>*- 000 


00 


CM 


m 


jt 


oc>o 


xo 


X 


0 


1 . 




U*“U 


•-U. 


CM 


m 


Jt 


c^ 


HU. 


u 


u 


i 




u u 


u 








■ 


uo 


u 


u 




45 



LL 











a: 




0 






C0 




• 






to 








— • 


* 














GO 






u 




to 




0 


< 




0 




• 


UL 




0 




♦ 


♦ 




CO 




0 






* 




z 


0 




3 C 






Z 




CO 




X 


N -4 




z> 




-J 


X 


CM 


to 




UJ 


u 


1 - 


0 




> 


UJ 


^0 


* 




* 


> 


^0 


3: 




• 


W 


V-QC 






ro 




< ♦ 


3: 




-f 


to— 


^GQ 


w 




• 


ocu 


CL< 


'S. 






h-<t 


3: 0 






W 


ZU- 


aDJcuJUJ 


CM 






->z 


Z>CCUX 


4 i 






♦ 0 


tOZJOO 


♦ 




• 


uu 


OiOZ * 


h- 




CM 


w l|( 


^o<z 


u 










< 




— H 


UJCO 




UL 




UJ 


UJtO 




♦ 




U 


CEO 


OXZ'^— 


CO 




NH 


-►u 




cM<in 




0 


— CD 


OCVIO!»CCMU. 


CD^O » 




z 


<*s. 


♦ * 3 CK-J 


tomuj 4 


• 


0 


K — 


<rr<in * ro< 


n<yy •» 


0 


0 


ujro 


1 cro 1 


-^0 •QUO 


CM C^ 




CD*- 


<<iLuinci: •<<QCCD^-^ 




• 


— 'to 


♦ 0 »— r- « 4 LLh- 


r- a 




OXCL 


OOXCM4'-''kOUJ< * to • 


•- z 


rOw( 


OOZ-' 


UJ QC r- r- X 0^ <l CO to bd < CM 


CM 






QN ♦ #winw 


--OUZ 1 


3: 


00 -f 


•4- to 


»<)< CM ro u. •- u. u. u •-» 0 UJ 0 0 0 


uaoc 


0 •« 


r-wO 


^b<C •wZ 4 l/)tO^XtOO • • • 


z • 


0^ 


II U.U 


00 <N-U-|-r-OO ♦ K ♦ »-»000 


►-‘OOO 


0 


><cz ♦ 


OOw II H-O •OOH-» 11 «U II II II 


3: II II II 


• 11 II bc:<K 


• II cxjcic<irn 11 


II xtocMZUJUJjTrn«juJUJ 4 rf) 0 »-HO 


KU 


*-|-<Oll IKOD II tOwOUUCM 


UJUUCM 


II to<o<<^ 


II uoto<i<totoh-<i II u*-^>-<i-o>*-<»-'i-ooa.u-cM 11 oc 



O 

OJ 

' fh 

in 

CM 



CD 

to 

a 

o 

CO 



to 

CL 

UL 

to 

a 

o 

« 



a: 

u- 



cc 

a 

in 

* 

K 

O 

< 

u. 

♦ 

o 

z 



ro 

in 

J 3 *<M 

h-CM 



z^uj II ikK0002:aiwa:-j«/^»-w jo:i/)H»-^-JZ »-• 

ZQC:Z»-»^UJUU<C7UJU-«^UJ»-iOU.UJ*-»‘-'Ot0UJOOt0w 

UJULOK<aO<GDt^VOCQi-<0>OOM>OOC!>aOOOQ. 



OCM 

^ m »CM 

in 3: •• 4 * •» 

CM -J •CM4 

^ UJ 4tOCMQ^ Jfc 

r— 

omtow in •^o 

i-ijTCLUL. *00 I UJO 

uoouto ^Z>OUJU • 
CLOtUOC UJI— troO<-^CM 

I •Q*- 03: I •--J I 

i-i I -fZ I-HQCIUOCUJUJ 

<>o 

CLVO’-' ♦0>0^0 

II Q.tOU< M M OC M HU 
i^ULCL II ULLUUJ^UJ4ZZ 
oo>-» ti ujzoooucMOa: 
C^QQOO*"* » '** ^*^h-OI3 

I (1.WI— Ot^tO*— ’H“0 Q 

' NHlktOXiOl— 

a»->OL> >o»-»oa*-»a:uiUJ 






CM ^ 

m 

4 



ou^ 

CM CM 



CMr- to 
CM CM 



46 




thesK66 

A computer program 



for forecasting the w 



3 2768 002 10687 4 

DUDLEY KNOX LIBRARY 













